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Abstract. We investigate qubit lasing in the strong coupling limit. The qubit is given by a 
Cooper-pair box, and population inversion is established by an additional third state, which can 
be addressed via quasiparticle tunneling. The coupling strength between oscillator and qubit 
is assumed to be much higher than the quasiparticle tunneling rate. We find that the photon 
number distribution is sub-Poissonian in this strong coupling limit. 



The coupling of an oscillator to a superconducting qubit [HE]; to a single-electron transistor 
[31 m [5], and directly to a tunnel junction [6] are promising possibilities to study quantum effects 
in electromagnetic or mechanical oscillators [3 [8]. On a principle level, all of those schemes can 
be used in multiple ways, for example for qubit read-out using the oscillator as detector [9], or 
the qubit can be seen as an artificial atom that can heat and cool the oscillator [10^ \TT\ I12j . 
This has been realized in two systems. A flux qubit has been used to to heat and cool an 
LC-oscillator [1'6\ [T^ , while lasing has been realized using a Cooper-pair box and quasiparticle 
tunneling to address a third state [TJl US] to create population inversion. We will investigate 
the single-qubit laser in the limit of strong coupling [T7]. This regime has an interesting lasing 
state with sub-Poissonian photon statistics. Our analysis extends previous studies in the strong 
coupling regime of the micro maser |18[ \W\ to a system where the atom is permanently coupled 
strongly to the cavity. 

We assume that our qubit is given by a Cooper-pair box coupled via the island charge to the 
oscillator. The qubit eigenstates are a superposition of the charge states = 0) and = 2) 
and are given by 

I t) = cos ^\N = 2)- sm^\N = 0) , | |) = sin ^|iV = 2) + cos ^\N = 0) . (1) 

The rotation angle is given by tan cf) = Ej /4:EcSNg, where Ej is the Josephson energy of the 
junction, Ec is the charging energy of the island, and 6Ng is the gate charge. The coherent time 
evolution of the system is determined by the extended Jaynes-Cummings Hamiltonian {h = 1) 

H = —AEaz + gicTz cos (f> + ax sin cp) (a^ + a) + uia^a . (2) 



A third state = 1) is involved in the lasing cycle, but here it is only included on the level 
of the master equation. In the strong coupling limit the time evolution of the density matrix is 
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Figure 1 . The average photon number (n) , the Fano factor F as a function of the frequency 
CO and the gate charge 5Ng- We see maxima in the photon number for the resonance condition 
AE = muj, where m = 1, 2, . . .. At resonance, the Fano factor has a minimum and the current 
has a maximum. The parameters are Ej/Ec = 0.18, A/Eq = 2.2, eV/Ec = 7, g/Ec = 0.01, 
K/{Ec/e^RT) = 0.1. 



given in sufficient approximation by a simple balance equation 

j 

where pi = {i\p\i), is the probabiHty of the system to be in the state rj_j=r^^ + Tj^f^ is the 
transition rate from state \j) to state \i) as obtained from the Golden rule, 

TZ = \mj)fi{E,, + eV) , r^^ = ^Y3^S^KW>^ (4) 

where Eji = Ej — Ei is the energy difference between initial and final state. We assume 
linear coupling of the oscillator to reservoir, x = + a, with coupling strength k/uj. The 
operator T = |1)(2| + |0)(1| decreases the charge of the island by one. The current trough a 
superconducting junction with resistance Rt is given by 

I{E) = di^fi^) (1 - /(^ + E)) N{u)N{lo + E) , (5) 

where N{E) is the dimensionless superconducting density of states and f{oo) is the Fermi function 
at temperature T. 

The crucial ingredient for any laser is the creation of population inversion. For the specific 
example we discuss here population inversion is created by quasiparticle tunneling. Without the 
coupling to the oscillator the quasiparticle tunneling rates are given by F^"^ ~ ^ sin^((/)/2), 
and F'l^P^ w F^^ oc cos^((/>/2) . For 5Ng > we get cos^{(t)/2) > sin2(</)/2). This means for the 
rates that Fi^-j- > F^-^i and T^ i > Fi |. Therefore the system is most likely to be in the state | f)- 
This creates population inversion in our qubit and can be used to generate lasing if the system 
is coupled to an oscillator. The energy for this process is provided by the transport voltage eV. 

We can numerically calculate the eigenstates as a combination of the Fock states of the 
oscillator \n) and the eigenstates of the Cooper-pair box = 1),| t), | \,). Equation ([3]) 
can be solved numerically in the stationary case pi = 0. We are interested in the average 
excitation of the oscillator (n) = (a^a) and the width of the distribution around this average, 
which defines the Fano factor F = ((n^) — (n)^)/(n). Results for T = can be seen in fig. 
([1]). As expected we observe maximal excitation if the oscillator is at resonance with the 
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Figure 2. We compare analytical and numerical solutions for the average photon number (n) 
and the Fano factor F. Black line: analytical solution, red-dashed: numerical solution. a)(n) as 
a function of k. The photon number decreases like k~^. For this plot we have chosen SNq = 0.1. 
b) The Fano factor F as a function of 5Ng. We see that the Fano factor becomes smaller than 
one for large 5Ng- For this plot we have chosen K/{Ec/e^R) = 0.1. For both plots we used the 
parameters Ej/Ec = 0.18, l^/Ec = 2.2, eV/Ec = 7, g/Ec = 0.01. 



energy difference AE. The Fano factor at these positions becomes especially small, which 
shows us that the density distribution has a sharp peak around the average value of the 
photon number (n). Directly at resonance the Fano factor can even be smaller than one, 
which means that we have a sub-Poissonian distribution. One can also observe higher order 
resonances, AE = mu, for m = 1,2,.... Close to this resonances our master equation is valid 

/ \m-l 

for g I (jry {n)/uj] ^ k, V/eRx, where m is the corresponding order of the resonance. Only if 

this condition is fulfilled the off-diagonal matrix elements of the density matrix can be neglected. 

To understand better the emergence of sub-Poissonian photon statistics we analyze the 
system at the resonances. On resonance the eigenstates are given by the dressed states 
|ib,77.) = (I t)|n)ib| l)\n+m))/\/^. We have now two sets of states, the even states |=t,n) and the 
odd states |l,n) = \N = l)\n). Quasiparticle tunneling leads to transitions between even and 
odd charges, this means it will cause transitions between the two sets. Oscillator dissipation 
does not change the charge of a state and therefore only causes transitions within each set. 
We can plug the eigenstates into the expression for the transition rates (jH). The quasiparticle 
tunneling rates between the dressed states are given by the same rates as tunneling between the 
eigenstates of the qubit as long as Eij <C eV. Assuming that (n) is large we can approximate 
the oscillator dissipation rates as 

r(i/±,n),{l/±,n-l) = in+ l)nK , rf^;±,nUl/±,n+l) = n (n + 1)/^ , (6) 

where n = {^e^/^sT _ Using theses rates we can write the balance equation Q as 

/0±,n = ^r'j'P^pi,„ + ^r'l^j^^/oi^n+m + (n + l)(n + 1)kp±,„+i + /innp±,„_i (7) 
- Qr;jP + ^r^P +{n + l)nK + n(n + 1)k^ /5±,„ , 

+(n + 1)k ((n + l)pi,n+i - npi^n) + nK (npi,„_i - (n + l)pi,n) . 

We can now find a solution for the average oscillator excitation (n). To do this we multiply 
eqs. ([7]) by n and sum over all n. We know that for sufficiently low dissipation the distribution 
is peaked around (n) S> 1, therefore we can neglect pi^. In the stationary limit we get a set of 



linear equations which can be solved for rij = ^n''^Pi,n ({^ ^ +> ~) !})■ The average number of 
photons is then given in good approximation by (n) = ?i+ + ?i_ + rii. This yields 



/pqp pqp _ pqp -pqp x 



+ n 



(8) 



We see that at resonance the oscillator excitations are inversely proportional to the dissipation 
rate k and linear in the order of the resonance m. The equation we derived here is rather 
general and describes any strongly coupled three- level laser. We can now use the same method 
to calculate higher moments of the photon number n. If we multiply equations (jT]) with 
and sum over all n it is straight forward to derive an approximate equation for (n^). For large 
photon numbers (F?''^ , T'j'']^ ^> F?''^, F^j^*!) the Fano factor reduces to 



Here we see that at the first order resonance (m = 0), low temperatures n ~ and for strong 
coupling, the Fano factor can be smaller than one. In fig. [2^) we compare the analytical results 
with the numerics. We have chosen 6Ng and k such that (n) is large enough to fulfill all 
approximations we have made to derive equations ([8]) and ([9]). In this case the numerical and 
analytical results fit perfectly. 
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